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1. Introduction 

The heavy-ion experiments at RHIC and soon LHC are probing deeper and deeper into the 
high temperature region of QCD. As such more and more experimental results are presented that 
require a more complete knowledge of finite temperature QCD. Naturally one turns to lattice QCD 
to give some answers to the questions arising. However some questions of interest are notoriously 
difficult to answer from a lattice perspective. This is in particular the case for the analysis of 
spectral properties of hadron correlation functions. Although here the main interest is in extracting 
information on hadronic spectral functions at low frequencies, this cannot be achieved without 
good control over the large frequency region. Conelation functions thus need to be controlled at 
short as well as large distances. High accuracy data on large lattices at several values of the cut- 
off are needed to control finite volume and lattice cut-off effects. Only then a reliable continuum 
extrapolation becomes possible which is a prerequisite for extracting dependable physics results. 

In this combined talk and poster contribution we present results from a lattice analysis of 
meson correlation functions at finite temperature [|[]. We emphasize here our systematic analysis 
of the finite volume and cut-off dependence of thermal meson correlation functions. This analysis 
leads to the conclusion that we indeed can take the continuum limit for the vector and pseudo scalar 
current correlation functions in a large Euclidean time interval. We also calculate several thermal 
moments of hadron spectral functions to better explore the spectral properties at low frequencies. 
We close with a discussion of the constraints arising from our analysis for the determination of the 
electrical conductivity in the quark gluon plasma at finite temperature [Q, Bh. 

2. Meson Correlation Functions 

Our key observable of interest is the Euclidean correlation function of a given particle current, 

J v = q(z,x)y v q(T,x) , (2.1) 

where choosing the appropriate gamma matrix we obtain the vector particle channels for y v = 
where pi = 0, ...,4 and the pseudo scalar for y v = 7s- 

In this work we analyze the Euclidean temporal two-point correlation function for the above 
channels at fixed momentum, 

Gp V (T,p)=Y t G lxv (z,x)e i P-*, (2.2) 

where 

G„ v {x,x) = {J fl (x,x)4(0,0)), (2.3) 

while we denote thermal expectation values with (...). In this work we will set the momentum to 
zero, p = 0, and suppress the second argument in the hadronic correlation functions. 

At high temperatures the spectrum of meson resonances is more and more changed by thermal 
effects; the width of resonances will broaden and the onset of the continuum in the spectrum may 
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shift. As a result the interest in the analysis of the above current-current correlation function shifts 
to their representation in terms of a spectral function p^, v 

r ur\ r d(0 n t n t\ cpsh(w(T-l/2r)) 

G ^ T) = .L 2^ P ^ (0) ' r) sinh(co/2r) • (2 " 4) 

From here on we denote with /iv = n the sum over the space-like components of the vector spectral 
function pa, while the full vector spectral function is given by py = Poo + Pa- We also use the 
notation PS for /iv = 55. Note that the correlation function Goo(T,r) is connected to the net 
number of quarks (q — q) in a given flavor channel, i.e., no(z) = J d xJq(x,x). As the net quark 
number is conserved, it does not depend on time, wo(t) = no. Thus the corresponding correlation 
function is constant in Euclidean time and the spectral function is simply given by, 

PooM = -2n% q (o8{(o) , (2.5) 

with the quark number susceptibility 

Xq = -j d 3 x(J Q (T,x)4(0,d)) . (2.6) 

Consequently the time-like correlation function obeys the relation Goo(t,T) = —% q T, which im- 
mediately leads to an exact relation between Gh(t, T) and Gy (t, T): 

G ii (T,T)=x q T + G v (z,T). (2.7) 

In the free field limit the spectral functions are seen to increase quadratically for large frequen- 
cies co. One obtains: 

3 

pa(ffl) = 2TiXaOo8(co) + — oo 2 tanh(oo/4T) (2.8) 

271 



3 

Pv{co) = poo(w) + pu{(0) = — w 2 tanh(w/4r) (2.9) 

2n 

3 

p PS (co) = — w 2 tanh(w/4r) (2.10) 
4n 

At finite temperature the exact cancellation of the 5-functions in the Pv{ooi) spectral function is only 
approximate, as poo(o) continues to be proportional to a 5-function due to the connection with net 
quark number, while the 5-function in Pu(oo) is subject to thermal effects. As a consequence this 
contribution is smeared out into a Breit-Wigner shaped peak [||] : 

pu((0) -> pl w ((0) = 2 Xq gw ffl2 ^ 2)2 +(! + «)• ^« 2 tanh(w/4r) , (2.11) 

where k parametrizes corrections to the high frequency behavior of the free field limit. This smear- 
ing in the low frequency region is directly related to the appearance of finite transport coefficients 
in the thermal medium, e.g. the electrical conductivity, 

G_ _ Cem Ph(®) 

T 6 £o^o coT ' v ' 

where C em is the sum over the squared charges of the contributing quark flavors. The pseudo scalar 
spectral function on the other hand does not per se contain an additional low frequency contribution; 
a connection with transport phenomena is not expected. 
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3. Thermal Moments of the Correlation Function 

A useful set of observables that helps to characterize the structure of spectral functions, is 
given by "thermal moments" of the spectral function at vanishing momentum. At a given order 
we define these moments as the Taylor expansion coefficients of the correlation function expanded 
around the midpoint of the Euclidean time interval 1 , i.e., around xT = 1/2, 

1 f°° dco /co\ 2 " Ph(co) 



(2n) 1 dG H (?T) 
Cr 



H ~ (2n)! d(tT) 2 " 



Tr=1/2 (2«)!7o 2n\TJ sinh(o)/2r) 



(3.1) 



and 



j \2n 



Gh (tT) = £g^ (--TT) . (3.2) 

n=0 V Z 



The infinite temperature, free field limit of the correlation function can be calculated analyti- 
cally. For massless quarks one obtains [||] : 

<«(,T) = 2<«(,T) = 6 T'U-2^-±^^ + 2^n\ (3. 3) 

V sm (2ktT) sin \2tzxT) J 

G{l ee (xT) = T 3 + G{ ree (xT) (3.4) 
Using these results the first three non- vanishing moments are then given by 



G (0),free = 2G (0),free = ^jree = 2ji ^ , . . ( 

r2 



G (2)Jree = 2G (2),free = 28^ ^ (3 6) 

G^ free = =2G$ free = l ^fT* . (3.7) 

In our lattice simulation we will also analyze ratios of moments. In the free field limit they are 
given by 

R^°) - ?V _ ^(2,0) _ 147^ _ 

K V,free ~ (0),free K PSJree ~ g - Z/ - OM 
try 

*&1 = 18123 (3-9) 

<jL = -^- >0.407, (3.10) 

which can be obtained from the ratio of the mid-point subtracted correlation functions and the 
corresponding free field values: 

MrT) _ Gh(tT) — GjP _ G% ] ( (4)2) d( 4, 2) u 1 



A£" (^) " (TT) - " V 1 + RHfrJ{ 2 + - 1 ' (3 - n) 

Note that the distinction between H = ii or H = V is unnecessary when evaluating Ah(tT), as dif- 
ferences in hte correlators, which arise from contributions of 5-functions in the spectral functions, 
are eliminated in subtracted correlators. 



1 As the spectral functions as well as the integration kernel in the spectral representation of the correlators are odd 
functions of the frequency, all odd moments vanish. 
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Table 1: Simulation parameters for the generation of gauge field configurations on lattices of size x N T 



4. Simulation Parameters and Results 

4.1 Computational Details 

Our numerical results are obtained from quenched QCD gauge field configurations generated 
with the standard SU(3) single plaquette Wilson gauge action Using an over-relaxed heat bath 
algorithm configurations were generated on lattices of size Ng x N x , where 32 < N a < 128 and 
N x = 16, 24, 32 and 48, with a separation of 500 updates per configuration. Correlation functions 
and plaquette expectation values calculated on these configurations have been checked and are seen 
to be uncorrelated. Details on the statistics collected are given in Tab. |. 

All calculations presented here have been performed at a single temperature value, T ~ 1.45r c . 
The gauge couplings have been adjusted accordingly for the different temporal lattice sizes. To do 
so we extrapolate known results for the critical coupling f5 c (N t ) and the square root of the string 
tension \fo This is achieved by fitting T/y/o using the ansatz suggested in []|] within the 
range of /3 G [5.6,6.5] and extrapolating to our relevant region of j8 E [6.8,7.8]. All simulation 
parameters are given in Tab. [j]. 

In the fermion sector we employ the clover improved Wilson action with non-perturbatively 
chosen clover coefficient csw [0] an d determine the correlation functions using an even-odd de- 



composed CG-algorithm. The hopping parameter K is chosen to be close to its critical value [11] 
and tuned to give approximately constant quark masses for four of the examined lattices. Here the 
quark masses are estimated using the axial Ward identity (AWI) to calculate the AWI current quark 
mass niAwi for the different values of cut-off. Whereby we use the non-perturbatively improved ax- 
ial vector current with coefficient ca noted in [|l2|]. In the next step we calculate the renormalization 
group invariant mass thrgi and rescale to the commonly quoted MS-scheme at the scale /I = 2GeV. 
The results are also listed in Tab. [T[ 

4.2 Numerical Results 

Due to the subtlety of the effects in question it is not very illuminating to examine the correlator 
directly. To illustrate this we show the correlation functions of the vector and pseudo scalar currents 
together with their analytically obtained free continuum and free discretized, lattice counterparts 
in Fig. |l|. Clearly the exponential decay of the correlation function dominates and obscures all 
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Figure 1: The vector (left) and pseudo scalar (right) correlation functions together with their free continuum 
(magenta lines) and free lattice (open symbols) counterparts. (Blow-ups): The ratio of free lattice and free 
continuum correlation functions. 



other interesting physics features. To circumvent this effect it is useful to look at certain ratios that 
largely cancel the exponential part of the Euclidean time dependence of the correlators. 

One of these ratios is that of the correlation function itself divided by the free correlation 



function given in Eq. 3.3, 



R H (r,T)= G f H{Tj) . (4.1) 
K ' G f H ree (r,T) 



The immediate advantage of a ratio as that given in Eq. |44] is that it gives a direct handle on the 
deviation of the correlation function calculated at finite and infinite temperature, respectively. As 
the continuum as well as lattice free correlation functions are known analytically, this also provides 
insight into the importance of lattice cut-off effects as function of Euclidean time. The insertions 
in Fig. [I] show the ratio of the free lattice and the free continuum correlation functions. They, of 
course, show that cut-off effects are largest at small distances. Moreover, it becomes obvious that 
with increasing cut-off values (larger A^ T ) the onset of cut-off effects is shifted to smaller distances. 



As we employ the local current correlations introduced in Eq. |2JJ all results have to be renor- 
malized multiplicatively, 



ff = {2KZ H )q(T,x)y v q(T,x), (4.2) 



where the renormalization constants have been determined non-perturbatively for the vector [[1 
and perturbatively up to two-loop order for the pseudo scalar channel [13]. They are listed in Tab. |]. 



In the case of the vector correlation functions the above statement holds true also for the time-like 
component Gqo{^,T) because the local lattice current is not conserved at non-zero lattice spac- 
ing. Consequently the time-like correlation function Goo(r,r) may not be strictly T-independent. 
However, as all vector current-current correlation functions are subject to the same renormalization 
constants a rescaling of the correlation functions Gy(%, T) and Ga(z, T) or even the ratio Rh{?, T) 
by Gqo(t, T) yields a quantity independent of renormalization. 
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Table 2: Quark number susceptibility (% q /T 2 ) calculated on lattices of size 128 3 x N z . The quark number 
susceptibilities have been renormalized using the renormalization constants listed in Tab. 1 . The last column 
gives the result of a continuum extrapolation taking into account cut-off errors of ff{a 2 ). 
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Figure 2: The time-like component of the vector spectral function, Goo('£T)/r 3 , calculated at T ~ 1.45r,;. 
The left hand part of the figure shows the volume dependence of Gqo{tT) /T 3 for N T = 16 and 32 < N a < 
128. The right hand figure shows the cut-off dependence of Gqq{xT) /T^ for Na — 128 and 16 < N z < 48. 



4.3 Finite Volume and Lattice Effects 

4.3.1 Quark Number Susceptibility 

As mentioned above the time-like component of the vector correlation function, Gqo(tT), 
will be used for rescaling in the following, so it is the first channel to be tested for finite volume 
and cut-off effects. We find these effects to be small, as can be seen from Fig. ||, where we plot 
—Goq(xT)/T 3 . The quark number susceptibility may then be extracted from the long distance 
behavior of the correlator, x q /T 2 = -G q(tT ~ 0.5)/Y 3 . 

On the left of Fig. || we show the results for the N r = 16 lattice with fixed cut-off while varying 
the spatial extent. Except for aspect ratio N /N t = 2 all results agree within statistical errors at the 
level of 1%. On the right the cut-off dependence at fixed spatial extent N a = 128 is shown. Here 
too systematic effects in the quark number susceptibility are seen to be small. 

Note even though we implement the non-conserved local current of the time-like vector cor- 
relation function, we obtain results that are T-independent to a very high degree. The deviations 
at small distances xT might be understood as lattice effects. The results for the quark number 
susceptibility are summarized in Tab. ||. 

4.3.2 Vector and Pseudo Scalar Correlation Functions 

Turning to the vector and pseudo scalar correlation functions we stress that the large range of 
available spatial lattice sizes at fixed cut-off, with aspect ratios ranging between 2 < N a /N x < 8, 
allows to quantify finite volume effects. Moreover, the large temporal extent of maximum N T = 48 
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Figure 3: The vector (left) and pseudo scalar (right) correlation functions, calculated on lattices of size 
N a x Nx at T ~ 1 A5T C , where "128 3 x 24 — 2" denotes the lower quark mass on this lattice. The results are 
normalized using the corresponding free correlation functions. Shown are data for tT > 1 /N z only. 



on an isotropic (!) lattice reduces the lattice spacing to O.Olfm at T ~ 1.45r f . The variation of 
N T by up to a factor of three allows to control lattice cut-off effects. In the following it will be 
shown that finite volume effects in the correlation functions are well under control and a controlled 
extrapolation to the continuum limit is indeed possible in a large Euclidean time interval. Note that 
that our current analysis improves on systematic errors that were present in earlier calculations of 
the vector spectral function performed by employing the same discretization scheme [14]. 

In Fig. |3| we show results for the vector and pseudo scalar correlation functions using the 
ratios introduced in Eq. 4. 1 . On the left hand side the ratio for H = V and on the right hand side 
for H = PS is shown for all available lattice sizes. Data sets with fixed spatial size at N a = 128 
while varying the cut-off N T shown in black. Data sets with fixed cut-off (N t = 16) and varying 
volume are shown in color. For one value of the cut-off (N r = 24) we performed calculations for 
two different values of the quark masses. We find that finite quark mass effects are small and well 
within 2%. From the fixed cut-off (colored) N t = 16 results in both plots finite volume effects for 
zT > 0.3 are seen to remain within a few percent even for the largest Euclidean time separation at 
zT = 0.5. As a consequence these results show that finite volume effects are under control. For 
zT < 0.3 in the left H = V plot it is immediately apparent that cut-off effects are large in the ratio. 

In the pseudo scalar case the situation concerning cut-off effects is not immediately evident. 
The ratio shown in Fig. ||(right) shows large deviations from the free field behavior even at short 
distances. At all distances the correlator thus seems to be controlled by large non-perturbative ef- 
fects. Moreover, the analysis of cut-off effects is obscured by the fact that data have been rescaled 
by using renormalization constants, which are known only perturbatively. This introduces unknown 
uncertainties. To eliminate at least these uncertainties we show in Fig. || the pseudo scalar correla- 
tion function normalized by the pseudo scalar correlation function at zT = 0.5. As we focus on the 
cut-off dependence we only show equal quark mass N a = 128 results. The left hand figure shows 
the pseudo scalar correlator normalized by free continuum correlator and in the right hand figure 
the free lattice correlation function has been used. From the left hand plot the cut-off dependence 
becomes immediately evident and we conclude that, similarly to the vector case, cut-off effects 
above zT = 0.3 are small and increase with decreasing zT. Even though their effect becomes ap- 
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Figure 4: The pseudo scalar correlation function normalized by the free continuum (left) and free lattice 
(right) correlation functions, both additionally rescaled by the correlator at xT = 1 /2, Gp$(xT = 0.5). 



parent for xT < 0.3, they do not dominate the behavior of the correlation function as in the vector 
channel. Actually the right hand side of Fig. ^ indicates that the T-dependence of the cut-off effects 
is similar to that of the free lattice correlation functions, as any cut-off effect is hardly visible also 
for xT < 0.3. 

4.4 Continuum Extrapolation 

The results shown in the previous section indicate that the continuum limit can be taken for the 
vector correlation functions at least for rr>0.25. To do so we use a quadratic ansatz in aT = 1 /N r 
to fit ratios of the vector correlator and the corresponding free field values. We normalize these 
ratios using the quark number susceptibility XqjT 1 and perform fits at fixed temporal extent xT . 
As cut-off effects are large on the N t = 16 lattice we will only use data from the N t = 24, 32, 48 
lattices, for reference we will however include the former in our figures. 

In Fig. | we show results of this extrapolation in l/N] for xT = 1/2, 7/16, 3/8, 5/16andl/4 
where we used the free continuum as well as the free lattice correlation functions for normalization. 
Note that we performed extrapolations for all Euclidean times xT available on the N t = 48 lattice. 
Wherever the smaller lattices fail to have a corresponding point in xT we interpolate using a spline 
construction. Subsequently the errors are then calculated using a Jackknife-method. In Fig. || this is 
the case for xT = 7/16 and 5 / 16. The figure reveals that the continuum limit can be cleanly taken 
and consistent results are obtained by using the free continuum and the free lattice normalizations, 
respectively. 

In the left hand part of Fig. ^ we show the results of the extrapolation in the vector channel 
as described above. Note that the largest deviation from the free correlation function occurs at 
xT = 1 /2. In fact, the established bending of the vector correlation function is crucial for a quanti- 
tative description of the low frequency region of the vector spectral function in terms of the ansatz 



suggested in Eq. |2.1 1| . The short distance part of the correlation function obviously can be well 
described by the free spectral function including the correction factor ( 1 + fc) as also has been done 
in Eq. [Ul| 

In the right hand part of Fig. ^ the corresponding result for the pseudo scalar correlation func- 
tion is shown. Here it is not possible to suppress the renormalization effects using suitable ratios 
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Figure 5: The ratio Ry(zT) normalized by the quark number susceptibility using the free continuum 
(full symbols) and lattice (open symbols) correlation functions over l/Nl Shown are the results of 
tT = 1/2,7/16,3/8,5/16 and 1/4, for legibility the results have been offset by 0.1 respectively. For 
TT = 7/16 and 5/16 spline interpolations were used to estimate the corresponding results on the N T = 24 
lattice. 




Figure 6: (left): The ratio Ry(rT) normalized by the quark number susceptibility using the free continuum 
correlation functions over xT and its corresponding continuum extrapolation, (right): The unnormalized 
Rps(tT) including its continuum extrapolation. In both cases the extrapolation was done as described, 
filling in spline interpolations when necessary. The solid curve in the left hand figure shows the fit discussed 
in section 5. 



of correlation functions. The extrapolation necessarily also includes this ambiguity. As the corre- 
lator normalized by its value at the midpoint was found to be almost cut-off independent and as 
finite volume effects were seen to be small renormalization effects dominate the uncertainty of the 
extrapolation. 

4.5 Curvature of the Vector and Pseudo Scalar Correlation Functions 

As discussed in section 3, thermal moments give additional insight into the spectral repre- 
sentation of hadronic correlation functions. They are especially interesting as they are obtained at 
the largest Euclidean time separation where the correlation functions are most sensitive to the low 
frequency region of the spectral function. In particular, the lower orders of the thermal moments 
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restrict the magnitude of the low frequency contribution to the spectral function and thus to the 
correlation function. 

In order to extract thermal moments we examine the quantity An{lT) defined in Eq. |3.1 ll 
Once more we rescale this ratio of subtracted correlators by the quark number susceptibility. In 
Fig. f7] we show results obtained from the vector and pseudo scalar correlators, respectively. We 
perform an extrapolation of Ah(xT) to the continuum, exactly as outlined for the correlation func- 
tions themselves. The extrapolated data is then fitted to a quartic polynomial as indicated in the 



Taylor-expansion in Eq. [3j] to obtain A#(tT) at xT = 1 /2. In the vector channel this gives 

G (2) 

,J* = 1.067 ±0.012 (4.3) 

r (2), free 

where H may denote either ii or V as noted in section 3. Combining these results with those of the 
continuum extrapolation for the vector correlation functions we obtain the ratios Ry 2 ® and R^'°\ 

R [ y fi] = 27.187 ±0.286 < R { yf ree and /?j?' 0) = 19.217 ±0.193 > R^f ree . (4.4) 

Repeating this analysis also in the pseudo scalar channel we obtain Aps(zT)/G^ shown in the 
right hand part of Fig. \j\ and the following results for the second moment, 

^3^(2) 

PS -0.7912±0.0012 .thus R& 0) = 10.932 ±0.017 < R { }?L„ . (4.5) 



°>S U PS 

These results reveal some interesting properties of the individual thermal moments; a combination 



of Eqs. |-k3j and |44j indicates that the second thermal moment is closer to the free field value for the 



H = V case and farther away for the H = ii, respectively. To evaluate Eq. 14^ in this way we need 



Gf s /Gf s free , which from Fig. § can be seen to be larger than 1. Subsequently, even though we 
are not able to extract the latter quantity without fully controlling renormalization effects, we can 
conclude the second moment must be closer to the free field limit than the zeroth. 

In both channels we also tried to examine the fourth thermal moment, but our numerical results 
unfortunately do not permit a conclusive determination of this value. 



5. Electrical Conductivity 

The results obtained for the vector correlation function and its continuum extrapolation, as well 
as the result on the second thermal moment put stringent bounds on the magnitude and shape of 
any contribution to the low frequency behavior of the vector spectral function. The small deviations 
from the free vector correlation function also suggest that the spectral function of the free theory 
is a good starting point for an analysis of the vector spectral function at finite temperature. We 
thus used as an ansatz for the spectral function the form given in Eq. [2.11| . This ansatz depends on 
four temperature dependent parameters: the quark number susceptibility % q {T), which we already 
extracted from the time-like component of the vector correlation function, the strength (cbw(T)) 
and width (T(r)) of the Breit-Wigner contribution and the higher order corrections to the high 
frequency free field spectral function, which we parametrize at present by a constant k(T). Already 
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Figure 7: (left): The mid-point subtracted vector correlation function normalized to the corresponding 
difference for the free vector correlation function. Shown is Av(tT) but normalized by the quark number 
susceptibility, (right): The pseudo scalar case Aps(tT) normalized by Gpg. The Fits in both Figures obey a 
quartic ansatz as indicated by the definition of A#(tT) and are shown within the interval tT 6 [0.2 : 0.5]. 



with this ansatz we obtain good fits for both the spatial (Go) and vector (Gy) correlation functions. 
In fact, a combined fit to the continuum extrapolated vector correlation function in the Euclidean 
time interval [0.25 : 0.5] and the second thermal moment, gives excellent results with a % 7 ' jdof 
below unity. For details on the fitting procedure and a more elaborate discussion of the results we 
refer to |JIJ]. The parameters obtained using this ansatz are: 

2c B wX q /r= 1.098 ±0.027 , Y/T = 2.235 ±0.075 , (1 + k) = 1.0465 ±0.003 . 

This fit is shown in Fig. |(left). Of course, as a consequence of this fit ansatz we also obtain a 
result for the behavior of the spectral function close to ft) = 0, i.e. we can deduce the electrical 
conductivity of the quark gluon plasma at T ~ 1.45T C : 

T 6 a>->o (oT 3 T v ; v 

We stress, however, that this result is a consequence of the particular ansatz used to fit the 
vector correlation function. An important question is, of course, to what extent this ansatz is unique 
or allows for modifications, in particular at low energies, which will influence the determination 
of the electrical conductivity. We intend to address this question by performing fits within a larger 
class of spectral functions as well as the Maximum Entropy Method [yj. 



6. Summary 

We have presented a detailed analysis of light meson correlation functions at T ~ 1.45r c in 
quenched QCD. For the vector current channels we find that finite volume and cut-off effects are 
under good control in a large Euclidean time interval. Here it is possible to take the continuum 
limit. The calculation of the second thermal moment and its inclusion in fits greatly helped to 
constrain the fit parameters. This led to an estimate of the electrical conductivity of the QGP at 
T ~ 1.45r o the significance of which requires further investigations in a larger parameter space. 
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In the pseudo scalar channel deviations from free field behavior are much more pronounced. 
In particular the analysis of cut-off effects is more difficult, as the perturbatively computed renor- 
malization constants introduce additional systematic uncertainties. However, rescaling the results 
by the pseudo scalar correlation function at the midpoint yields a largely cut-off independent re- 
sult. Also finite size effects are found to be small in the pseudo scalar channel. This suggests that 
a spectral analysis of the pseudo scalar correlation functions should yield reliable results for its 
frequency dependence and may suffer only somewhat from an imprecise knowledge of the overall 
normalization. We will address the spectral analysis of the pseudo scalar correlator elsewhere. 
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